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Abstract 

There are two paradigmatic frameworks for treating quantum systems coupled to a dissipative 
environment: the Caldeira-Leggett and the Ambegaokar-Eckern-Schon approaches. Here we recall 
the differences between them, and explain the consequences when each is applied to a zero dimen¬ 
sional spin (possessing an SU(2) symmetry) in a dissipative environment (a dissipative quantum 
dot near or beyond the Stoner instability point). 
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The diagrammatic technique for non-equilibrium systems developed in the pioneering 
works of Schwinger and Keldysh plays a predominant role in theoretical condensed mat¬ 
ter physics 12 . It is designed to tackle real time evolution of systems at and away from 
equilibrium. Following the developments of the last two decade^®, it now provides a non- 
perturbative tool to tackle interaction induced strong correlations in quantum many-body 
systems. In this paper we discuss an important prototypical problem, a quantum zero di¬ 
mensional degree of freedom in a dissipative environment, in which the Keldysh technique 
is of tremendous use, providing insight into the physics involved. 


I. GENERAL PERSPECTIVE 

We consider here the dynamics of a quantum system coupled to a dissipative environ¬ 
ment. The resulting equation-of-motion is stochastic, which can be formulated on any of the 
following three levels: (i) a fully classical Langevin equation, where both the variables are 
classical (expectation values of observables) and the frequency range of interest is Tioj < ksT. 
For Ohmic dissipation the noise spectrum is white; (ii) a semi-classical hybrid description, 
within which the variables are still classical coordinates, but one acknowledges the fact that 
the noise may be quantum, having high frequency component, Tiu> > ksT (Ref. 6 -^); (iii) a 
full-fledged quantum mechanical description, according to which the noise may contain high 
frequency quantum components, and the variables of the quantum Langevin equation are 
operators within the Heisenberg description. This approach is practiced, say, in the field of 
quantum optics'. 

A paradigmatic framework to present a dissipative environment, in a way that connects 
to our preformed classical intuition, is to model Ohmic resistor quantum mechanically. We 
mention here three approaches: 

1. The Caldeira-Leggett (CL) modelling 8 : One introduces an effective circuit consisting 
of an L-C transmission line (with infinitesimal imaginary term), that may extract energy 
and current from the bare quantum system, (cf. Fig. [1^) 

2. The Ambegaokar-Eckern-Schon (AES) modeling 9 10 : Here we model a tunnel junction 
(cf. Fig. [l)}) , assuming explicitly that its transparency is low, hence only lowest order 
contributions in the tunneling should be accounted for. The resulting Hamiltonian represents 
reservoir degrees of freedom that give rise to dissipation. Traditional applications of the CL 
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FIG. 1. Three approaches to envisioning dissipation: a) an LC circuit (wave-guide), extracting en¬ 
ergy from the system; b) a dissipative tunnel junction, suitable for the AES picture; c) the Landauer 
picture consisting of a tunnel barrier (or a tunnel juction) coupled to dissipative reservoirs. 


picture employed extended coordinates (this, however, is not a must; the CL action in the 
case of a spin degree-of-freedom consists of compact coordinates). By contrast the AES 
approach introduces compact (periodic) coordinates. 

3. The Landauer picture 11 13 . Here one models the resistor by a tunnel barrier (of 
arbitrary transparency) (cf. Fig. [TJ^ for the single channel case). The contribution of this 
tunnel barrier to the resistance is given by R/( 1 — R ), where the reflection probability off 
the barrier is equal to the modulus square of the reflection amplitude, R— |r| 2 . This elastic 
backscattering process yields the magnitude of the resistor; the actual inelastic dissipation 
takes place in the connected reservoirs. Such a model has been discussed, for example, in 
Ref. 11 . We shall not consider this picture here. 


The outline of this paper is the following: in Section[TT]we will briefly review earlier works, 
emphasizing the difference between the CL and AES approaches to dissipative dynamics, 
focusing on charge dynamics. The gauge symmetry underlying charge transport is U(l). In 


Section III we will recall the physics of a quantum dot (QD) tuned to be near (but below) 
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the Stoner instability. As such, the QD supports large magnetization. Ignoring fluctuations 
in the magnitude of the spin, the spin degree of freedom possesses an SU(2) symmetry. The 
coupling of such a QD to external leads gives rise to dissipation, which is formulated and 
studied within the framework of the AES (Section IV). In Section [V] we first compare our 
AES analysis for the spin case to our results obtained within the CL framework. We then 
note that this AES vs. CL contrast differs from the AES vs. CL in the standard charge 


U(l) case. We conclude in Section VI 


II. CALDEIRA-LEGGETT VERSUS AMBEGAOKAR-ECKERN-SCHON: THE 
CHARGE U(l) CASE 

We consider the dynamics associated with current through a resistor, and compare the 
two paradigmatic representation thereof: CL and AES. 


A. CL action 

The CL action of a current biased linear resistor (modeled as a transmission line) reads 

iS CL = - J dhdt 2 a(t h t 2 ) + source (!) 

Here the dimensionless phase variable ip(t) represents the effective flux variable <L(f) via 
tp = 27t($/* 1 * 0 ), where $0 — h/e is the flux quantum. The voltage across the resistor is given 
by V = and is the degree-of-freedom canonically conjugate to the charge that 

has flown through the resistor Q — f dtl. In (Jl]) a(ti,t 2 ) is the kernel of the Ohmic bathP. 
Dropping the time-local terms (important for avoiding renormalization of the non-dissipative 
part of the action) we obtain 

iS C L = J dtidt 2 a(ti, CiMHMO) + iS source . (2) 

Note that in Keldysh notation this action may be written as 



( 3 ) 
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where c,q refer to the classical and quantum components on the Keldysh contour J15 . The 
subscripts R, A, K refer to the retarded, advanced, and Keldysh components of the matrix. 

Employing the relation between the retarded and the advanced components of the kernel 
a, a A (t 2 ,ti) = — Q! (fi, f 2 ) we may write the action as 


iScl — iScL T iSch + *5 


with 


iSc L = 2i J dtidt 2 -t 2 )\ <p q (h)(p c (t 2 ) , 

iS CL = \ J dtidt 2 a K (ti - t2)(Pq{ti)<p q {t 2 ) , 


and 


lS source = l J dtIex ( t )j± <p q {t) , 

One maj! 6 rewrite the Keldysh term of the action, employing the decoupling 

e ti s cL — I e % f dtn t( t )<P g ( t ) e h I 


It follows that 


m)^)) = Ti a K (t 1: t 2 ). 


(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


The resulting Langevin equation-of-motion is obtained by calculating the variation i5ScL/S(p q (t) 
0. The equation obtained is 

Mt) 


R 


— Rx + 5I(t) , 


( 10 ) 


where SI(t) = e£(i) represents stochastic current noise. We note that the noise is additive , 
and is not affected by the bias current. 

In deriving Eq. 10 we have used the fact that the dissipative bath has an Ohmic spectrum! 8 , 
implying that 

11 h ' 2 _ 11 h 2 

( 11 ) 


R .. 1 1 h 2 „ R/ . 11 h 2 

Im a R (t) = - - 5 (t) or Re a R (u) = - - j :<u , 


where R is the resistance. The variation over the retarded part of the action leads to 


iSSnr 2vr . /' 

CL = — 2i / dU 


s$ q (t 1 ) <h 0 


Im a H (t 1 — t 2 ) 


2n- 


< hn 




( 12 ) 


The Fourier transform of the current noise correlator is given by 


{6I(t 1 )6I(t 2 )) u = -a K (u) 
h 


( 13 ) 
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At equilibrium 


a K (co) = a R (co) — a A {u) 


coth 


hu> 

2 knT 


(14) 


The fluctuation-dissipation theorem follows from Eqs. (12) and (14). 

.. T/ , .. hu , Tioj 

(6Ilti)6I{t 2 )) u = / , coth^-^ . 


(15) 


We note that the additivity of the noise and its independence of the bias current (Eq. 10) 
imply that the noise is independent of I ex , i.e., absence of shot noise. 


B. AES action 

The AES action now is given by 

iS AE s = ~ J dt-idt 2 a(ti,t 2 ) (1 - cos [<^(fi) - ip{t 2 )]) + iS source . (16) 

The source term is the same as in the previous case. Similarly to the CL case, Eq. Q, one 
may write the action as 


iSaes = + iSaes + is 


( 17 ) 


The retarded part is essentially identical to that in the CL case, having to do with the fact 
that t\ and t 2 are very close to each other (cf. Eq. @), which allows us to expand the 


cos(...) term in Eq. (16). The Keldysh term, though, is very different: 

iS AES = \\ dt\dt 2 a K (ti - t 2 ) 

{[cos(^(ti)] g [cosp(i 2 )] 9 + [sin(^(t 1 )] 9 [sin^(t 2 )] 5 } , 

Decoupling the action, employing two auxiliary fields, and £ 2 , one obtains 10 

e s S AES = J X)^X)^ 2e s/ dt?i (?ih)[cos^(d] 5 +€2(i)[sm¥i(t)]J 

,1 J dtidt2h [ aK ] t t (^1 (*1 )€l (*2)+C2 (*1 )^2(*2 )) 


xe‘ 


(18) 


(19) 


The resulting equation-of-mot ion for the AES action is 

nr = ~ e6 sil1 ( 2ir fj;) + cos 


( 20 ) 
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This equation can be cast into the form of Eq. (10) by writing dl(t) = dli(t) + dl 2 {t) with 
the two independent terms of current fluctuations defined as 


$ 

dli = — e£i sin ( - 

Vq 


$ 

dl 2 = e £ 2 cos ( 2tt —- 
Vq 


( 21 ) 


The equation-of-motion (20) implies that the noise is non-additive, as can be shown 
explicitly from the following iterative procedure. The zeroth iteration gives <f> c = Vt, where 
V = I ex R. Next we introduce a correction <f> c = Vt + 5$ c and obtain 


^c(t) 

R 


= —e£i sin 27T 


Vt + d$ c 


+ e £ 2 cos 27r 


Vt + d$ c 

$n 


( 22 ) 


The first iteration consists in dropping 5$ c in the r.li.s. of Eq. (|22j). The resulting stochastic 
terms give rise to shot noisd 10 -^ (unlike the CL equation-of-motion). For eV ksT we find 


(fill(ti)dli (t 2 ))cj—> o — (^L 2 (ti)h/ 2 (t 2 )) a ,_ s .o — x e 77 

Z K 


(23) 


III. A QUANTUM DOT NEAR THE STONER PHASE TRANSITION 

Over the past few decades the physics of quantum dots has become a focal point of 
research in nanoelectronics. The introduction of the Universal Hamiltonian 16 19 has made 
it possible to take into account the effects of electron-electron (e-e) interaction within a QD 
in a controlled way. This approach is applicable for a normal-metal QD when the Thouless 
energy Exh and the mean single particle level spacing 5 satisfy #qd = Exh/d » 1. Here 
$qd is the dimensionless conductance of the QD. The single particle level spacing is given 
by d r\_/ l/(Uz/ 0 ), where V is the volume of the QD and u 0 is its density of states (DoS) 
and therefore d ~ 1 /L d for a d-dimensional QD. The Thouless energy, Exh, is the inverse 
time-of-flight (or diffusion time) of an electron across the quantum dot. 

Within this scheme interactions are split into a sum of three spatially independent contri¬ 
butions in the charging, spin-exchange, and Cooper channels. Ignoring the latter (see below) 
the charging term leads to the phenomenon of Coulomb blockade, while the spin-exchange 
term can drive the system towards the Stoner instability 20 . In bulk systems the exchange 
interaction competes with the kinetic energy leading to Stoner instability. In finite size 
systems mesoscopic Stoner regime may be a precursor of bulk thermodynamic Stoner in- 
stabilitjJ™: a new phase, intermediate between paramagnetic and ferromagnetic emerges, 
in which the total spin of the QD is finite but not extensive (i.e., not proportional to the 
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volume of the dot). The mesoscopic Stoner regime can be realized in QDs made of materials 
close to the thermodynamic Stoner instability. 

A quantum dot in the metallic regime, $qd 1, is described by the universal Hamilto- 
iiiarf® 

H = H 0 + H c + Hj + H\. (24) 

The noninteracting part of the universal Hamiltonian reads 

H 0 = J2 , (25) 

a,a 

where e a denotes the energy of a spin-degenerate (index cr) single particle level a. The 
charging interaction term 

H c = E c (N - N 0 f (26) 

accounts for the Coulomb blockade. Here, Ec = e 2 /(2 C) denotes the charging energy of 
the quantum dot with the self-capacitance C, N 0 represents the background charge, and 
N = J2a,a a a<r a aa i s the operator of the total number of electrons of the dot. For the 
isolated quantum dot the total number of electrons is fixed and, therefore, the charging 
interaction term can be omitted. The term 

Hj = -JS 2 (27) 

represents the ferromagnetic (J > 0) exchange interaction within the dot where S = 
J2 a a a a 1 S<Jia 2 a a <y 2 ^ ie °P era t° r of the total spin of the dot. Here S CT1(T2 = (l/2)a (Tl(T2 , 
where a = (a x , cr y , cr z ) is a vector made of Pauli matrices. The interaction in the Cooper 
channel is described by 

Hy = \T'T, T=J«„ Al . (28) 

a. 

In what follows we do not take into account H\ for the following reasons. For the dots 
defined in 2D electron gas the interaction in the Cooper channel is typically repulsive and, 
therefore, renormalizes to zero 1 '. In the case of 3D quantum dots realized as small metallic 
grains, the interaction in the Cooper channel can be attractive, giving rise to interesting 
competition between superconductivity and ferromagnetism 21 23 . In that case we assume 
that there is a weak magnetic field which suppresses the Cooper channel. 


The starting point of our analysis of a dissipative Stoner QD (near the Stoner instability 
point) accounts for the QD Hamiltonian 

H dot = J2 e a a l,a a a,* ~ JS 2 , (29) 

a,a 

In doing so we ignore possible correlations between the charging state and the spin config¬ 
uration of the QD 24 . 

We note that for isotropic spin exchange interaction (Heisenberg model) the mesoscopic 
Stoner phase extends over 1/2 < J/5 < 1. For the anisotropic cas d 25 * 26 ! the lower boundary 
of this inequality slides towards 1, with no mesoscopic Stoner phase for Ising spin 24 27 . For 
the isotropic case the ground state spin S is the integer value (for even number of electrons 
on the QD) or half-integer value (for odd number) that is closest to <7/2(5 — J). This value 
increases with increasing J and diverges for J —)• 5, which marks the onset of the macroscopic 
Stoner ferromagnetic phase. Seemingly the problem is easy to tackle theoretically. The 
interaction terms of the universal Hamiltonian consist only of zero mode (zero wave-number) 
contributions, which commute with each other. The inclusion of the exchange term renders 
the problem non-trivial though: the resulting action, which consists of Pauli matrices, is 
non-Abelian (more specifically, it is underlined by an SU(2) symmetry). Attempts to study 
the problem from different points of view included the Ising limit 24 , perturbation theory in 
the Ising anisotropy 25 . An exact solution that employs states classified by the total number 
of electrons and the total spin 28 30 requires the calculation of Clebsch-Gordan coefficients 
which is not an easy task. In this way Alhassid and Rupp have found an exact solution for 
the partition function in the absence of Zeeman splitting. Elements of their analysis were 
then incorporated into a master equation analysis of electric 28 29 and thermal properties. 
Independently, a study of electron transport through a QD for low temperatures (T 5) was 
made in referencd- 32 . That analysis, accounting for the charging and exchange interactions, 
employed a master equation approach as well. 

An exact solution based on the Wei-Norman-Kolokolov approach had been presented in 
referencd^, and was then extended to include randomness-induced spectral fluctuations 34 . 
The tunneling density of states and the spin susceptibility were calculated; other thermody¬ 
namic and linear response correlations are calculable as well. The study of shot noise near 
the Stoner point was reported in 35 . 

We note that the exact solution approaches mentioned above, while elegant and power- 
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ful, are very difficult to generalize to more complex setups, in particular, to setups where 
external leads are added - a common mean for the introduction of dissipation. An efficient 
approximation, which can be generalized to such setups, employs adiabatic approximation 
of the spin stochastic dynamic^. 

IV. AES APPROACH FOR SU(2) SPIN 

Our approach 3 ^ can be viewed as a generalization of the Landau-Lifschitz-Gilbert (LLG)- 
Langevin equation 38 ^, central to the field of spintronics 4 -^, to a regime where quantum 
dynamics dominates. Stochastic LLG equations have been derived in numerous publications 
for both a localized spin in an electronic environment (a situation of the Caldeira-Leggett 
typef 1 42 and for a magnetization formed by itinerant electron^- 43 ®. In all these works the 
precession frequency was assumed to be lower than the temperature or the voltage, thus 
justifying the semi-classical treatment of the problem. In this regime the geometric phase 
did not influence the Langevin terms. 

Our derivation here is technically close to that of Ref. 43 . However, in contrast to Ref. 4 ®, 
we do not limit ourselves to small deviations of the spin from the instantaneous direction, 
but rather consider the action on global trajectories covering the entire Bloch sphere. 

To demonstrate the emergence of an AES-like effective action we consider a quantum 
dot with strong exchange interaction coupled to a normal lead. The Hamiltonian reads 
H = H dot + Hi ead + H tun . The quantum dot is described by the magnetic part! 4 ^ of the 
universal Ilainiltoniari 13 


H do t = e a a l,a a a,a ~ J S ' 2 + B S , (30) 

a,cr 

where S = (1/2 ) a a,<r 1 <T <r 1 ,<y 2 a aa 2 4S the operator of the total spin on the quan¬ 

tum dot, B is the external magnetic field, and J > 0 is the corresponding “zero mode” 
ferromagnetic exchange constant. The Hamiltonian of the lead and that describing the 
tunneling between the dot and the lead are standard: Hi ead = X) 7 ,<r e 7 c 7 ,<r c 7 a and H tun = 
Sa, 7 ,cr h/, 7 a a,o- c 7 >o-‘h^‘ c " We assume here a non-magnetic lead. Here 7 is the orbital quantum 
number describing eigenmodes of the lead. 

We consider the Keldysh generating functional Z = J HTUT exp [i iSvi/, where the 
Keldysh action is given by bR = §K dt (?'T<9 f T — H) (plus the necessary source terms which 
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are not explicitly written). Here, for brevity, T denotes all fermionic fields and the time f 
runs along the Keldysh contour. After standard Hubbard-Stratonovich manipulation^' 33 36 
decoupling the interaction term — JS 2 we obtain Z — f DA 4 exp [i Sm] , and the action for 
the bosonic vector A d(t) reads 


iSm = tr In 


r - 1 

' Jr dot 


-v 


-0 G 


-i 

lead 



M | 2 

4J 


(31) 


are 


Here G dot = [id t - e a - ( M(t ) + S) • cr/2], while G Iead = - e 7 . Both G0 ot and G lead 

matrices with time, spin, and orbital indexes. We introduce M(t ) = A't(f) + B. Expanding 


(31) in powers of the tunneling matrix V and re-summing we easily obtain 
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iSm = tr In 


n- 1 
^lead 


+ tr In 


r <-i 

^ dot 


dt 


\M — B\ 


4 J 


(32) 


K 

where the self energy reads £ = VG{ ead V^. The first term is trivial, i.e., it would never 
contain the source fields. Thus, it will be dropped in what follows. 

Rotating frame. We introduce a unit length vector 


n(t) = (sin 9 cos (j), sin 9 sin 0, cos 9) (33) 

through M(t ) = M(t)n(t ) and transform to a coordinate system in which n coincides with 
the z-axis n(t) ■ cr = R(t)a z R'(t). This condition identifies the unitary rotation matrix R 
as an element of SU(2)/U(1). Indeed, if we employ the Euler angle representation 


R = exp [—(*0/2)0%;] exp [-(i9/2)a y \ exp [-(iif/2)a z \ , (34) 


then the angles 0(f) and 9{t ) determine the direction of n(f), while 0(f) is arbitrary, i.e., the 
condition n(t) • <r = Ra z Rf is achieved with any value of 0(f). Thus, 0 represents the gauge 
freedom of the problem. We introduce, first, a shifted gauge field y(f) = 0(f) + 0(f). This 
way a periodic boundary condition, e.g., in the Matsubara representation R(r) = R(t + 0), 
is satisfied for \{r + /3) = y(r) +47rm (The fact that m is integer is intimately related to the 
spin quantization 16 ). We can always assume trivial boundary conditions for y, i.e., m = 0. 
We keep this representation of the rotation matrix R also for the Keldysh technique. 

We perform a transition to the rotating frame and obtain 


iS M 


= tr In 


R] ( G 


-i 

dot 




K 


BM 

2 J 


(35) 
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(we omit the constant term oc |.B| 2 ). For the Green’s function of the dot this gives R)G d ° t R = 
idt — e a — M(t ) cr z /2 — Q, where we define the gauge (Berry) term as Q = R^(—id t )R = Q\\ + 


6 a y — (j) sin 9 a x 


Q. l. Here Q\\ = [0(1 - cos 6) - x] a z /2 and Q± = - exp [?'yoy] 

Note, that Q depends on the choice of the gauge held x■ Finally, we obtain 

'M 2 BM 


iSm = tr In 


G~ d L -Q- -R'sfl 


— I 


dt 


K 


4J 


2 J 


exp [i(j)a z }/ 2. 


(36) 


where G^ tz = id t -e a - (1/2) Af(t) cr 2 . 

To hnd the semi-classical trajectories of the magnetization we need to consider paths 
M(t), 6{t ), 0(t) on the Keldysh contour such that the quantum components are small. The 
quantum (g) and classical (c) components of the fields are expressed in terms of the forward 
( u ) and backward ( d ) component^, i.e., 4> q (t) = 4> u (t) — 4>d(t) and 0 c (t) = (4> u (t) +4 > dif))l 2- 
Performing the standard Keldysh rotation 4 we thus obtain 


iS m = tr In 


+ i dt 


r 

^ dot,z 


Q - R ] tR 


2 J 


— i dt 


M c M q 
2 J 


(37) 


where G~^ tz = T x G~^ tz . The local in time matrix helds Q[t) and R{t) also acquire the 
2x2 matrix structure in the Keldysh space, e.g., Q = Q c t x + Q q To/2, where t x ^ i2i0 are the 
standard Pauli matrices. 


The adiabatic limit. Thus far we have made no approximations. The action (37) governs 
both the dynamics of the magnetization amplitude M(t) and of the magnetization direction 
n(t). Here we focus on the case of a large amplitude M (more precisely, M fluctuates 
around a large average value M 0 . Such a situation arises either on the ferromagnetic side of 
the Stoner transition or on the paramagnetic side, but very close to the transition. In the 
latter case, as was shown in Refs. 3 ^®, it is the integration out of the fast angular motion 
of n which creates an effective potential for M, forcing it to acquire a finite average value. 
More precisely the angular motion with frequencies oj max [T, B] (we adopt the units 
h = kg = 1) can be integrated out, renormalizing the effective potential for the slow part 
of M(t). The very interesting question of the dissipative dynamics of slow longitudinal 
fluctuations of M(t ) in the mesoscopic Stoner regime will be addressed elsewhere. Here we 


focus on the slow angular motion and substitute M{t ) = M 0 . Thus, the last term of (37) 
can be dropped. We note that in the adiabatic limit we may neglect Q± as it contributes 
only in the second order in dn/dt 36 . 
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The idea now is to expand the action (37) in both Q (which is small due to the slowness 
of n(f)) and WT,R (which is small due to the smallness of the tunneling amplitudes). A 
straightforward analysis reveals that a naive expansion to the lowest order in both violates 
the gauge invariance with respect to the choice of y(f). One can show that the expansion in 
WT,R is gauge invariant only if all orders of Q are taken into account, that is if (Gdot Z ~Q )~ 1 
is used as zeroth order Green’s function in the expansion. This problem necessitates a clever 
choice of gauge, such that z ~ Q) -1 is as close as possible to Gdot,z , he., the effect of Q 
is “minimized”. 


Choice of gauge. As the action (37) is gauge invariant we are allowed to choose the most 


convenient form of y(f). We make the following choice 

Xc{t) = (1 - cos 60(f)) , 

Xq{t) = <i> q {t){l- cos 9 c (t)) , 

which satisfies the necessary boundary conditions, i.e., y 9 (f = ±oo) = 0. 


(38) 


Here we present a detailed justification of the gauge which is presented in Eq. (38). Ideally 


we should have chosen a gauge that would lead to Q\\ = 0. Seemingly, this might have 
been achieved with the choice y(f) = (j>(t) (1 — cos 9(t)) on both branches of the Keldysh 
contour. This choice, however, violates our desired boundary conditions as the integrals 
over x accumulated between t = —oo and t = +oo on the upper and on the lower Keldysh 
branches are different. Such a difference would show up as non-trivial boundary conditions 
on x q at either f = —oo or t — +oo. In other words, had we selected y(f) = 0(f) (1 — cos 9(t )) 
we should have violated the requirement y 9 (f = ±oo) = 0. We note, though, that to linear 
order in the quantum components the condition y(f) = 0(f) (1 — cos 0(f)) yields x q = 0 9 (1 — 

t f • . -I 

cos 9 C ) + 9 q sin 9 C 0 C , leading to y g (f) = / dt' <f> q (t')(l — cos0 c (f')) + 9 q {t') sin0 c (f') 0 c (t') = 


0?(f) (1 cos 0 c (f)) + f dt' sin 9 c (t') 9 q (t') 0 c (f') - 9 c (t') (j> q (t') 


The first term vanishes at f = 


±oo but not the last term. We thus include only the first term in y g , leading to Eq. (38). The 


gauge (38) satisfies the boundary conditions and leads to the desired cancellation = 0, 


whereas the quantum component of Qy remains nonzero: 


Q\\,q = g ° z Sill6l c ~ 0c(fq 


(39) 


At the same time this choice allows for the expansion of the Keldysh action in the small 0 ? 


and 6 q as there are no 0 g terms remaining in (39) 
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Berry phase (Wess Zumino Novikov Witten (WZNW) action). Expanding the zeroth 


order in £ term of the action (37) to first order in Q we obtain the well known in spin 
physics (see, e.g., Refs.® 4 ') Berry phase (WZNW) action 


iSwzNW — — n J dt tr 

which after a straightforward calculation reads 


Gdot.zMQ iuW 


iS 


WZNW 


= iS clt sin 9 r 


4*c@q d c (j)q 


(40) 


(41) 


where S = N(M 0 )/2 is the (dimensionless) spin of the dot. Here N(Mq) is the number 
of orbital levels of the dot in the energy interval M 0 around the Fermi energy. Roughly 
S = M 0 p dot /2 , where p dot is the density of states averaged over the energy interval M 0 . The 
effects of mesoscopic fluctuations of the density of states were considered in 


AES action. The central result of the current paper is the AES-like 9 effective action, 
which we obtain by expanding (37) to the first order in R)T,R: iSAES = —tr Gdot,zR t £ R 
This gives 



(*i—* 2 ) 


Rc(t 2 ) 

Rq{t 2 ) 
2 


(42) 


where g = (h /e 2 )(G7f- + G 4 ,)/ 2 . Here G a = 2tt (e 2 /%) \V\ 2 p lead p a dot is the tunneling conduction 
of the spin projection a, p'J^ are the densities of states at the respective t and j. Fermi levels, 
whereas the density of states in the lead, pi ea d, is spin independent. The standard ) 10 Ohmic 
kernel functions are given by a R (u>) —a A (u>) = 2u and a K (u> ) = 2u coth(a;/2T). The action 
(42) strongly resembles the AES action 1 -®, with U( 1) exponents exp [f(p/ 2 ] replaced by the 
SU( 2) matrices R. Fixing the gauge of R is an essential part of our procedure. 


Semi-classical equations of motion. From the effective action (42) we derive the semi- 
classical equation of motion. We follow the ideas proposed in Ref. 6 . Using the representation 
R = A 0 a 0 + iA x a x + iA y a y + iA z a z , with A 0 = cos 


e 

2 

cos 

X 

2 

, A x = sin 


sin 

XI C4 

1 

© 


A y = — sin 


cos 


, A z = — cos 


Sill 


| we rewrite the AES action (Eq. (42)) as 


SaES = SaeS + S AESi where 

iSaes = - 2 ig J dtidt 2 lma R (ti - t 2 j\ J2 A j,g(G)A jjC (t 2 ) , 


( 43 ) 
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and 


iS AES = -| / dtidt 2 a K (ti - t 2 ) A i,q(h) A i,q(h) 


( 44 ) 


Here j = 0,x,y,z. The Keldysh part of the action (44) leads to random Langevin forces. 


This can be shown !- 6 using the Hubbard-Stratonovich transformation 


iS K 

£ LC> AES = 


n oi, 

\j=0,x,y,z 


X 


exp 


dt\i Zj A m \ + iS Z 

j=0,x,y,z 


(45) 


where the action S £ is given by 


J dtidt 2 

“ i 


a 


K 


-i 


(ti-ta) 




(46) 


In other words, = 8jkg& K (ti — t 2 ) and (4) = 0. We obtain the Langevin 


equations Eq. (J47]) from 5iS to tai /S(j) q (t ) = 5iS to tai/ dO q (t) = 0 , where iS tota i = LSb-MSm/zw + 
iS% ES +j dt Yj i^j.\j, r Here iSs = —iS'y B J dt sin 9 C 6 q is the action related to the magnetic 
held (in ^-direction). Prior to performing the variation of the action, the held \ is replaced 


according to the gauge fixing choice (Eq. (38)). Finally, we use a R (t) = (d t + C)5(t) (the 
constant C is important for causality but drops in our calculation) and obtain the following 
equations of motion: 


d c + g sin 6 c (j) c = t] 6 , 


sin 6 r 


7 B) - gO c = r) (t> . 


(47) 


Here g = ^ and 7 = ( JPdot ) 1 is the “gyro-magnetic” constant of order unity. The Langevin 
forces (torques) are given by 

1 


Ve = 




2S COS 2 


2S Sm 4 


2S COS i 


1 


04 

2 

0c 

2 

0_ c 

2 

0 C 


4 COS 0, 


Xc 


+ iy sin 


Xc 


2 S' Sm 2 


c Xc . c . Xc 
iz COS — + 4o S111 — 


4 sin ( 0 C - y ) - 4 cos 


t • Xc , Xc 
4 Sill y - 4 cos y 


Xc 


(48) 
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The l.h.s. of Eqs. (47) represent the standard Landau-Lifshitz-Gilbert (LLG) equation^ 38 " 
(without a random torque). The r.h.s. represent the random Langevin torque. The latter is 
expressed in terms of four independent stochastic variables (j = 0 ,x,y,z), which satisfy 
= 5jkga K {ti — £ 2 ) and (£j) = 0. On the gaussian level, i.e., if fluctuations of 


9 C and <f c are neglected in Eqs. (48), the Langevin forces go and are independent of each 
other and have the same autocorrelation functions: ( 770 (^ 1 ) 770 (^ 2 )) = 0 and ( 770 (£ 1 ) 770 (£ 2 )) = 
( 77 ^(£ 1 ) 77 ^,(£ 2 )). We emphasize that, in general, the noise depends on the angles 9 C and 


(j) c leading to complicated dynamics within Eqs. (47). In the classical domain, i.e., for 
frequencies much lower than T, we can approximate (Ci(G)^fc(O)) = 4gT<5 (£ 1 — £ 2 ) 5jk • Then 
= ( 770 (£ 1 ) 770 (£ 2 )) = (gT/S 2 )5(ti — £ 2 ). Thus, the situation is simple and we 
reproduce Ref. 39 . 

Effective temperature. In the quantum high-frequency domain the situation is different. 
We cannot interpret the four independent fields as representing the components of a 


fluctuating magnetic field. A close inspection of equations (47) shows that in the regime of 
weak dissipation, S 1 and g • Cl, the spin can precess with frequency B = at an 
almost constant 9 for a long time of order (shorter than) (gB)^ 1 . For such time scales we can 
approximate </> c = Bt and Xc = (1 —cos 9 c )<f> c = (1 —cos 9 c )Bt. Thus the Langevin fields f n in 


(48) are multiplied by fast oscillating cosines and sines with frequencies oj cos = B cos 2 (6 ) c /2) 
and to S i n = B sin 2 (0 c /2). Thu^ 8 


(li,Ati) r l 4 ’At 2 ))u=o = -^[cos 2 (0 c /2)ajj Aa,) + sin 2 (9 c /2)ajf (w. 


(49) 


In the quantum regime T B these correlation functions differ substantially from the 
classical ones, (q^itfri^it')) ^=0 = {ve{t)ve(t'))u}=o = gT/S 2 . Thus, if the spin could be held 
for a long time on a constant 9 C = 9 0 trajectory (one possible way to do so was proposed in 
Ref.® 1 ), the diffusion would be determined by the quantum noise at frequencies co c and co s , 
which are governed by the geometric phase. More precisely, the spread of 9 C and 0 C (in the 
rotating frame) will be given by (A 9) 2 = sin 2 9q (A f) 2 = Dt , where 

D = (g/S 2 )T efI , (50) 


and the effective temperature is calculated from (49) to be 


T eff = \ cos4 (y) coth 


B 

— Sill 
2 


4 M coth 


B_ 2 9^ 
2T C ° & y 2 / 

J0o\ 

: Sill — 


B 
2T 


(51) 
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We emphasize once again that this semi-classical analysis is valid for a highly non-equilibrium 
situation is which the spin is driven and is kept artificially at a trajectory with 9 C = 9 0 ^ 0. 
Semi-classical approximation. We are now ready to discuss the physical meaning of the 


semi-classical approximation, i.e., the expansion of the action (42) up to the second order 
in 9 q and <f > q . The non-expanded action is periodic in both 9 q and <j ) q . The periodicity in (f) q 
corresponds to the quantization of the z spin component S z = S cos 9 C . By expanding we 
restrict ourselves to the long time limit, in which S z has already ’’jumped” many times by 
A S z = 1 in the course of spin diffusion. We neglect, thus, higher than the second cumulants 
of spin noise (see, e.g., Ref. 49 for similar discussion of charge noise). We obtain, however, a 
correct second cumulant with down-converted quantum noise (similar to shot noise in the 
charge sector). This is due to the ’’multiplicative noise” character of our Keldysh action 
(42) similar to the original AES casd^ (see alscf 50 () . 


Equilibrium dynamics near 9 C = 0. In the absence of external driving at T R, Eqs. (47) 
lead to fast relaxation of the spin towards the north pole of the Bloch sphere, i.e., 9 C = 0. 
Here we show that the effective temperature introduced above looses its meaning in this case. 
Near the north pole the spherical coordinates are not adequate and we rewrite the Langevin 
equations (47) in cartesian coordinates. Namely, we define x = sin 0 C cos </> c ~ 9 c cos(J) c and 
y ~ 9 c sh\(f) c . The new Langevin equations for x and y (valid for i,i/< 1) read 


x = -By - gBx + 2S ^ + ~ 2 ^ (£* “ 9tv) , 

V = Bx — gBy + ^ (£„ + gQ • (52) 

A straightforward analysis of these linear equations leads to the stationary widths (standard 
deviations) of order Ax = Ay ~ 1 /\/S. Taking into account the standard relation ( S 2 ) = 
(Si) + (S^) + (Si) = S(S + 1), we observe that in the pure state S z = S the following relation 
holds (Si) + (S^) = S 2 (Ax 2 + Ay 2 ) = S. Thus, fluctuations of order Ax = Ay ~ l/VS 
are purely quantum (they would be of this order also for AS Z ~ 1) and the semiclassical 
analysis is inapplicable in this case. 


V. CL VS. AES 


In this Section we compare the SU(2) AES model described in Section IV with the 
straightforward generalisation of the Caldeira-Legget model for the spin SU(2) case. We 
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further notice the similarity between {the difference between AES and CL in the U(l) case} 
and {the difference between AES and CL in the SU(2) case}. 


A. CL in the spin SU(2) case 


The Caldeira-Leggett action arises from the interaction of the type H int = h • n. Here 
n = S/S' and the vector field h represents isotropic fluctuations of an effective magnetic field 
with the Keldysh correlation function (' TKh n (t\)h m (t 2 )) = ga/H,^) Ai,m, where the times t\ 
and t 2 are on the Keldysh contour. The hied h can, in reality, be due to, e.g., the Kondo 
coupling of the localized spin S to the electron-hole continuum. The coupling constant g 
is chosen so that the equations of motion are exactly the same as in the AES case, where 
g was proportional to the tunneling conductance. Assuming the fluctuations are Gaussian 
one obtains the following effective action 


Scl = f / dti <f dt 2 a(t 1 ,t 2 ) (1 - n(H)n(t 2 )) 
2 Jk Jk 


(53) 


The Keldysh analysis similar to that presented above produces again equations (47), however 
the Langevin terms look different: 


Ve 

9<t> 


— (-^sin0 + 4cos0) , 

sin 9 cos 9 , . 

~2S~ ^ 2 _ ~2S~ ^ x C ° S ^ + ^ SU1 ^ 


(54) 


Only three random fields (n = x, y, z ) are needed. Their fluctuations are given by 
{€n(ti)£ m (t 2 )) = 8 nm goi K (ti — t 2 ). Exactly these equations are derived in Ref. 39 before 
making the high temperature approximation, which makes the cos (f> and sin (j) factors in the 
right hand side unimportant. In contrast to the AES case we obtain 


(Ve(ti)ve(t2))oj =o = 
(?M(n)?fa(o))„= o = 


9 

4:S 2 

9 

AS 2 


a K (u = B) , 
cos 2 9 C a K (u = B^j 


+ sin 2 6 C a K (u 



(55) 


We observe that the diffusion is not isotropic is this case. That is, in the ^-direction the 
diffusion is characterized by Dq = ( g/S 2 )Tg , where Tg 
we obtain D$ = ( g/S 2 )T/ with | cos 2 9 B coth 
We observe that Tg > T^. This anisotropy is most pronounced for 6 = 7t/ 2 and T B. 


= A B coth 


B_ 
2 T 


B_ 

2T 


For the ^-direction 
+ sin 2 9 T = cos 2 9 Tg + sin 2 9 T. 
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Once again we emphasise that the above mentioned diffusion takes place in a highly non¬ 
equilibrium case of a spin being artificially held on a trajectory with constant 0/0. At 
equilibrium, as in the AES case, the semi-classical analysis is not applicable. 


B. Comparisons: CL and AES for U(l) and SU(2) cases 


Here we compare the similarities and differences between the CL and the AES pictures 
in the U(l) charge case, with those in the SU(2) spin case. In the U(l) case the semi- 


classical equation of motion can be cast in the form of Eq. (10) both for the CL and the 


AES models. Yet, the Langevin term, i.e., the fluctuating current, SI, is entirely different 
in the two models at low temperatures, k&T <C eV. In the CL case 81 = e£ is produced by 
one stochastic variable £, whose noise spectrum is Ohmic at equilibrium. In the AES case 


one needs two independent stochastic variables £i and £2 (see Eq. 21). Both these variables 


have equilibrium Ohmic noise, yet, due to the multiplicative oscillating factors in Eq. (21), 


the noise of 51 at zero frequency is determined by the noise of £ 1j2 at frequency V. This 
leads to the appearance of shot noise. 

Analogously, in the SU(2) spin case, both CL and AES models lead to the semi-classical 


stochastic LLG equations of the form (47). The two Langevin terms (spin torques) rjg and 


are, however, different in the two models. In the CL case rjo and 77 ^ can be expressed (see 


Eq. 54) via three independent stochastic variables £ x , C y , £ 2 (all having Ohmic equilibrium 


noise spectra). In comparison, in the AES case one needs four independent stochastic 


variables £ x , £ y , £ z , £ 0 with Ohmic equilibrium spectrum (see Eq. 48) 


In both SU(2) CL and SU(2) AES models the noise is multiplicative. That is, in both 


Eq. (54) and Eq. (48), the independent Langevin variables are multiplied by trigonometric 


functions of the Euler angles 6 and </>. Thus, in both models, the frequency shifts are similar 
to those leading to the shot noise in the U(l) case. Yet, these frequency shifts are very 
different between the CL and the AES cases. We consider again the example of the spin 
being held artificially at the trajectory with 9 = 9 0 7 ^ 0 and precessing with frequency 
B. In the CL model the spectrum of £ 2 is not shifted, whereas the spectra of £ x and £ y 
are shifted by B. By contrast, in the AES case, the spectra of £ x and £ y are shifted by 
u cos = Hcos 2 (6 , o/2) and the spectra of £ 2 and £0 are shifted by u> S i n = B sin 2 (# 0 /2). Both 
these factors are geometrical and are determined by the Berry phase of the spin’s trajectory. 
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VI. CONCLUSIONS 


In this paper we review the Caldera-Leggett (CL) and the Ambegaokar-Eckern-Schon 
(AES) approaches to dissipation. We first remind the reader about the well known physics 
of dissipative charge dynamics underscored by the U(l) symmetry. Then, we provide an 
analogous treatment for the dissipative SU(2) spin dynamics. In both cases the Keldysh 
technique allows deriving the quasi-classical Langevin equations of motion. Except in the 
CL U(l) case, the noise is multiplicative, which leads to the admixture of the high frequency 
(quantum) noise components to the low frequency dynamics. This gives rise to shot noise 
in the charge dynamics and well as to the novel phenomenon of geometric dephasing in 
dynamics of large spins. 
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